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Abstract 

Analytical solutions of variable coefficient nonlinear Schrodinger equations 
having four-dimensional symmetry groups which are in fact the next closest 
to the integrable ones occurring only when the Lie symmetry group is 
five-dimensional are obtained using two different tools. The first tool is 
to use one dimensional subgroups of the full symmetry group to generate 
solutions from those of the reduced ODEs (Ordinary Differential Equations), 
namely group invariant solutions. The other is by truncation in their Painleve 
expansions. 

1 Introduction 

The purpose of this paper is to classify solutions of a general class of variable coefficient 
nonlinear Schrodinger equations (VCNLS) of the form 

i^Jt + f{x,t)tp^^ + g{x,t)\ij\^i: + h{x,t)i) = 0, 
/ = /i + ^/2, 9 = 9i + i92, h = hi+ih2, (1.1) 
f„g„h,eR, J = 1,2, /i^O, gi^O. 

with the property that they are invariant under four- dimensional Lie symmetry 
groups. This class of equations models various nonlinear phenomena, for instance 
see [1] and the references therein. Symmetry classes of (11.11) are obtained in [2] and 
canonical equations admitting Lie symmetry algebras L of dimension 1 < dimL < 5 
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are presented there. A suitable basis for the maximal algebra L (dimL = 5) is 
spanned by 

T = dt, P = d^, W = d^, B = td^ + ^xd^, D = tdt + ^xd,-^pdp, (1.2) 

which is isomorphic to the one- dimensional extended Galilei similitude algebra gs(l). 
Here ip & C is expressed in terms of the modulus and the phase of the wave function 

=p(x,t)e*'^(^'*). (1.3) 

An equation of class fll.ip admits this algebra as long as the coefficients /, g and h 
can be mapped into 

/ = 1, g = e + ig2, e = ±1, 5(3 = const., h = (1.4) 

by point transformations. This is nothing but the standard cubic nonlinear 
Schrodinger equation (NLSE). For the form of the coefficients obeying the constraints 
imposed by the Painleve test, we had been able to transform fll.ip to the usual NLSE. 
In two recent papers |3] and the conditions imposed by the Painleve test were 
shown to be equivalent to those having a Lax pair. Therefore these conditions are 
also necessary for integrability. 

We intend to present a detailed analysis of solutions to VCNLS equations in 
the absence of integrability using two different approaches. We focus on the 
canonical equations which are representatives of the equations from class fll.ip having 
four- dimensional Lie algebras. A list of four-dimensional symmetry algebras and the 
corresponding coefficients for the canonical equations is given in Table [1] (We hereby 
correct an error on L5 in f2\ and slightly improve the form of the equation). 



Table 1: Four dimensional symmetry algebras and the coefficients in (11.11) 



No 


Algebra 


f 


9 


h 


Conditions 


Li 
L2 
L3 

L5 


{T,D,,C,,W} 
{T,P,B,W} 
{T,P,D2,W} 
{P,B,D2,W} 
{P,B,C2,W} 


1 
1 

1 + ^/2 
1 

1 


e + i'j 
e + i-y 
e + i'y 

e+i-y 


(hi + ih2)^ 
ih2 


' t 

2hi+i{2h2-t) 


h2^0 
/2 7^0 
h2^0 




2(l+t2) 





7, hi, /i2 are constants and e = ±1. Wave function written in the polar form (II. 3p . 
the basis elements for the symmetry algebras are given by 

T^g„ P = 9., w^g., B = + (1.5) 

and 

11 1 

Ci = t^dt + xtd^ - -tpdp + -x^d^, Di = 2tdt + xd^ - -pdp 
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for Li, D2 = tdt + \xdx — \pdp for L3 and L4, and C2 = {\ + t^)dt + xtdx + \x^ for the 
algebra L5. We note that Li is a non-solvable and L2 is a nilpotent algebra and the 
other three are solvable and non-nilpotent. In addition, Li and L3 are decomposable 
whereas the others are not. 

These canonical equations do not pass the Painleve test for PDEs, therefore they 
are not integrable and will be the main subject of this study. We are going to apply two 
different methods: Symmetry reduction and Painleve truncated expansions. The first 
is to make use of the one-dimensional subalgebras of the four- dimensional algebras 
given in Table [1] and the second is to find a valid truncated series solution to the 
equation. 

The paper is organized as follows. In Section 2 we find the group-invariant 
equations for the canonical equations having four dimensional symmetry algebras. 
Section 3 is devoted to the analysis of the reduced systems and completes the study 
of the invariant solutions. In Section 4 we apply the method of truncated Painleve 
expansions to the canonical equations to obtain exact solutions. 

2 One-dimensional subalgebras and reductions to 
ODEs 

As we are interested in group invariant solutions we only need one-dimensional 
subalgebras. This is the case because we restrict ourselves to subgroups of the 
symmetry group having generic orbits of codimension 3 in the space {x, t} x {p,uj}. 

The classification of one-dimensional subalgebras under the action of the group 
of inner automorphisms of the four-dimensinal symmetry groups is a standard one. 
We do not provide the calculations leading to the conjugacy inequivalent list of 
subalgebras. The classification method can be found in for example in O El [7]. 

The main result is that every one-dimensional subalgebra of the symmetry algebra 
is precisely conjugate to one of the subalgebras given in the Table [2l 

Using these subalgebras, we perform the reductions leading to the ODEs. We 
exclude the subalgebras whenever the invertibility requirement is violated. This is 
only case for the gauge symmetry W and it does not lead to any group-invariant 
solutions. We first write the wave function in the form (11. 3p and obtain (11. ip as a 
system of two real second order nonlinear PDEs, given by 

~ pujt + fi{pxx- P^l) - f2{'^Px^x + P^xx) + 9iP^ + hip = 0, (2.1a) 

Pt + f2{Pxx - P(^l) + h{'^Px(^x + P(^xx) + g2P^ + h2P = 0. (2.1b) 

In this system coefficient functions with indices 1,2 are real and imaginary parts 
of /, g and h. They are all functions of x and t. For example, if we would like to see 
the system for the algebra Li, looking at the Table [1] we simply replace hi{x,t) of 
(12. lap by ^ and for the algebra L5 by jt^s, this time hi being a constant. 
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Table 2: One-dimensional subalgebras of 4- dimensional algebras under the adjoint 
action of the full symmetry group 



Algebra 






Subalgebra 


a,b,c e M, ei = =f1 

7 7 — 7 J. 1 


Li 


Li.i 


= {T + Ci + aW} 


L1.2 = {Di + bW} 


L^,3 = {T + cW} 


L2 


L2.1 


= {P} 


L2.2 = {T + aW} 


L2.3 = {B + bT} 




L2A 


= {W} 






L3 


L3.1 


= {T} 


L3.2 = {P} 


L3.3 = {T + eiW} 




L3A 


= {P + e,W} 


L3.5 = {D2 + aW} 


L3.6 = {T + eiP + bW} 




L3.7 


= {W} 






L, 


L4.1 


= {P} 


L4.2 = {B} 


L4.3 = {P + ei5} 




Laa 


= {D2 + aW} 


U.^ = {W} 






L5.1 


= {P} 


L,.2 = {B} 


L,.3 = {C2 + aW} 



Invariant surface condition for a specific subalgebra gives the similarity variable 
for the functions p and u. Use of this variable in (12. ip therefore reduces the number 
of independent variables in the system from two to one, converting it to a system 
of ODEs. These nonlinear systems of ODEs arise as first or second order nonlinear 
equations. 

In the following, each time we encounter a first order equation, we provide its 
solution right away. The task is harder for the second order equations, because we 
have to decouple the reduced system of equations before treating them. For specific 
values of the constants appearing in the reduced equations, we were able to succeed 
in decoupling the systems and left the search for solutions to the Section 3. 

2.1 Non-Solvable Algebra Li = {T, Di, Ci, W} 

Commutators for the basis elements of the four-dimensional algebra Li satisfy 

[T, D,] = 2T, [T, C] = Di, [Di, Ci] = 2Ci (2.2) 

with W being the center element, that is, commuting with all the other elements. 
The algebra has the direct sum structure 

Li = sl(2,R)©{Vr}. (2.3) 

Representative equation of the algebra is 

#i + ^xx + (e + + {hi + ih2)\^ = (2.4) 

X x^ 
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with the real constants e = =i=l, 7, /ii, /;,2. In polar variables it takes the form of a 
system 

e h 

- put + pxx- P^l + -p^ + = 0, (2.5a) 

7 s ^2 

Pt + 2pxOJx + poJxx + -p + -ijP = 0. (2.5b) 
2.1.1 Subalgebra Li.i = {T + Ci + aW} 

Invariance under the subalgebra Li.i implies that the solution has the form 



ip{x,t) = — ^exp 



t (a arctan t + ^7^^ + ^(0 



4(1 + 



^2 



We substitute p{x,t) = ^^^^^ and uj{x,t) = a arctan t + ^^j^^*^^ + P(^) in (12. 5p and 
obtain the reduced system of equations satisfied by M(^) and P{^): 



eM^ + [-^- -a^ + hijM- A^'^MP'^ + A^'^M" = 0, (2.7a) 

+ h2M + SC^M'P' + A^'^MP" = 0. (2.7b) 

We first need to decouple these equations to solve for the functions M and P. If 
(I2.7bp is multiplied by M and written as 

7M^ + h2M^ + 4^2(M2p')' = 0, (2.8) 

it is seen that an integral of (12. 8 p can be obtained for three different cases of the 
constants. 

(i) The case 7 = 0, /;.2 7^ 0. 

We define Y = Y{^) in ([23D as 

(M^py = = r, (2.9) 

which then easily integrates to 

M^P' = Y + C, C = const. (2.10) 

so that 

h2 



= -^^^Y'. (2.11) 



Hence the expression 

h 2{Y + C) 
4^2y 
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obtained in terms of Y has to be substituted into fl2.7ap . But the other terms in fl2.7ap 
must also be expressed in terms of Y, which is seen to be possible if one considers 
(I2.1ip . As a result we obtain a third order nonlinear ordinary differential equation 
from fl2TiD : 



Yr'''--Y'''+-Y'r'+—i^-^-a^+h)Y''--r'-^^{Y+cy = o. (2.13) 

From a solution of this equation, one can find the functions M, P of fl2.6p as 

m) = 2^/^, ^(0 = /^^^- (2-14) 

(ii) The case 7 7^ 0, /;,2 = 0. 

This time the function Y is defined in fl2.8p as 



{M'PJ = = r (2.15) 

leading to the integral 



M^P' = Y + C, C = const. (2.16) 

so that we have 

M^ = ^^Y'. (2.17) 

7 

From f l2.7ap the following third order ODE for Y is obtained: 



YY'" - ^-Y'" + ^-Y'Y" - — -p + r + ^ y'^ + ^(y + c)^r' + ^F'^ -- = 0. 
4 4 ^ 4 4/ 4 V 7 

(2.18) 

For a solution 1^(4^) of this equation M and P are going to be evaluated from 



M(o = ( - -ey) , pio = 11 (2-19) 

(iii) The case 7 = /12 = 0. 

In this case we can easily decouple the reduced system of equations. Integration 
of (ESD gives 

M^P' = C, P(0 = f ^d^^ C = const. (2.20) 



and from fl2.7ap we obtain the equation for M 

M" = -^M' + + - /ii)m + C'M~\ (2.21) 
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2.1.2 Subalgebra L1.2 = {Di + bW} 

Group-invariant solutions for subalgebra L1.2 will have the form 



ip{x,t) = exp 



i(61nx + P(0) 



X 

T 



(2.22) 



It is straightforward to see that M{C,) and P{C,) must satisfy 

eM^ + (^^-b^ + hi^M + ^(^ - Ab)MP' - A^^MP'^ + A^^M" = 0, (2.23a) 

7M^ + (/i2 - 26)M + ^(46 - ^M' + 8^'^M'P' + A^^MP" = 0. (2.23b) 
If we multiply fl2.23ap by M, we can write it in the form 

26 1 



7M^ + h2M' + C 



M^i — -- + AP' 
^ 2 



(2.24) 



Integration of fl2.24p is possible in three different cases, 
(i) The case 7 = 0, /i2 7^ 0. 

The function F(^) is defined in fl2.24p as 



We integrate to find 
and with the equality 



M^(--- + AP' 
2 

1 h 

P' 



e 

h2Y + C 



8 2^ 4^2 Y' 
h2 



(2.25) 
(2.26) 
(2.27) 



fl2.23ap can be completely expressed in terms of Y 



1. 



Y'Y'" — -Y"'^ + -Y'Y" + 
2 e : 



2^^,^^„ l/3 + 4/ii 2b 1 



i 4 



^-]Y' 



rl2 



2/lr 



r^-^(r + C)2 = 0. (2.2^ 



If a solution F of this equation is known, then P and M can be obtained from fl2.26p 
and dMZD- 

(ii) The case 7 7^ 0, /;,2 = 0. 

In (EMD the function Y{Cj is defined by 



.2/26 1 



+ AP' 

V ^ 2 



= >^ • 



We integrate fl2.29p equation to obtain 

1 b Y + C 
8 2^ 4^ 

7 



(2.29) 



(2.30) 



and hence 

M^ = --^^Y' (2.31] 
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of which substitution in fl2.23ap gives rise to an equation in terms of Y: 



rr"-5r«+iKY"+(l^-l+l)y-+^(K+c)V'+lry-^ = o. 

(2.32) 

P and M are going to be found from fl2.30p and fl2.3ip . 
(iii) The case 7 = /i2 = 0. 

If (I2.24P is integrated once, we get 

and substitution of this resuh in fl2.23ap leaves us with the second order equation for 
M: 

M" = C^M-^ - -^M^ - -i-(3 + 4/ii - 2b^ + ]^^)M. (2.34) 
44 lo^^^ 4 

2.1.3 Subalgebra L1.3 = {T + cW} 

Invariance under the subalgebra L1.3 implies that the solution will have the form 

il){x,t) = M{x)exp i{ct + P{x)) (2.35) 

and here M(^), P{$,) satisfy the system 

exM^ + {hi - cx^)M - x^MP'^ + x^M" = 0, (2.36a) 
7xM3 ^ }^^^ ^ 2x'^M'P' + x^MP" = 0. (2.36b) 

Similarly, can be arranged as 

-fxM^ + hiM"^ + x'^ {M^P')' = (2.37) 

and with arguments similar to the preceding algebras we obtain the following results, 
(i) The case 7 = 0, /i2 7^ 0. 

M{x) and P{x) are found from 

M{x) = (^-IxV^ ^ , P{x) = -h2 J^^^dx. (2.38) 
Here Y{x) satisfies a third order equation 



Y'Y'" - -Y"^ + -Y'Y" + 2('^ - - _ _^ ^ Q_ (2.39) 

2 X \x^ ) h2 x^ ^ ^ ^ ^ 



(ii) The case 7 ^ 0, /i2 = 0. 

M{x) and P{x) are to be evaluated from 



M{x) = ^ ' ^(^) = y (^ + ^)\/~^'^^- ^2-^°) 

is a solution to the equation 

Y'Y'" - + ^y/y// + ( ^^^^ " ^ _ 4 \ y/2 + l2(y + c^)2y/ + ^y'\I--Y' = 0. 
4 2x V 4x2 / X X Y 7 

(2.41) 

(iii) The case ^ — h2 — 0. 

M{x) is the solution of the equation 

M" = C^M-^ + f c - ^ W - -M3 (2.42) 

and P{x) is going to be evaluated from 

P{x)^J^dx. (2.43) 

2.2 Nilpotent algebra L2 = {T, P, P, M^} 

Nonzero commutation relation is [P, S] = ^W. The algebra contains the three 
dimensional abelian ideal {T, P, 1^}. The action of B on this ideal can be represented 
by the nilpotent matrix N 

/ 1/2 

\ 

In this case the canonical equation has the form 

#t + i^xx + (e + ^7)1^1^^ + ih2ilJ = (2.44) 
with the real constants e = ^1, /i2 7^ and 7. 

2.2.1 Subalgebra L2.1 = {P} 

The group-invariant solution of L2.1 has the form 

ijjix, t) = M(t) exp {iP(t)) (2.45) 

and M, P must satisfy 

eM^ - P' = 0, (2.46a) 

7M3 + /i2M + M' = 0. (2.46b) 







[T,P] 


] I T 


[W, B] J 


' \ W 
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We immediately integrate these equations and find 



M{t) = ( Miexp(2/i2t) - 



^1/2 



P{t) 



27 



In (^Mi - gexp{-2h2t) 



with arbitrary constants Mi, Pi. 



Pu 7 7^ 
7 = 



2.2.2 Subalgebra L2.2 = {T + aW} 

A solution invariant under the algebra L2.2 must be in the following form 

^/j{x,t) = M{x)exp i[at + P{x)) 
Here M, P have to satisfy 

eM^ -aM - MP'^ + M" = 0, 

+ h2M + 2M'P' + MP" = 0. 



For 7 = we multiply fl2.49bl) by M to obtain 

h^M^ + {M^py = 

and define the function Y{x) in this equation as 

{M^P'Y = -h2M^ = ¥'. 

Substitution of 

/Y + C 
y, dX 



in f l2.49al) leads to decoupling of the reduced system 

-Y''-' - 2aY''^ - 2hl{Y + = 0. 



Y'Y'" — -Y""^ — —^'3 

2 h2 



P is determined by fl2.52p and M is given by 



M{x) 



h2 
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2.2.3 Subalgebra L2.3 = {B + bT} 

(i) The case 6 7^ 0. An invariant solution of L2.3 is obtained in the form 

t2 



Functions M and P are solutions to the system 



^M^ - -^M - MP'^ + M" = 0, 
jM'^ + /12M + 2b'^M'P' + b'^MP" = 0. 



We can arrange (I2.56bl) as 

7M^ + h2M^ + b\M^P'y = 
and for 7 = define Y{^) such that 



2 D'V 



{M'P 



52 



Im^ = Y' 



from which we get 



Hence we have decoupled f l2.56al) in the form 



1 9f r 2/)2 

Y'Y'" - -Y"^ - ^Y"> - ^Y'^ - -^{Y + Cf = 0. 

2 /12 



P is obtained from (I2.59P and M is given by the formula 

62 



M(0 



/l9 



.y/\l/2_ 



fii) The case 6 = 0. 



= M(t)exp ^(- + P(t) 



(2.55) 

(2.56a) 
(2.56b) 

(2.57) 

(2.58) 
(2.59) 

(2.60) 

(2.61) 
(2.62) 



is the form of the group-invariant solution and the reduced system of equations is 



eM^ - P' = 0, 



+ (/i2 + ^)M + M' = 0. 
We solve this system by standard methods and get 

— 1 /2 

M(t) = f texp(2/i2t)(Mi + 27 / ^^P(-2^2t) ^^^^ 



(2.63a) 
(2.63b) 



(2.64) 
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and 

P{t) = \ ^ ^ ' (2.65) 

I IvIr/^^^^^^^ + ^i' 7 = 0. 

In the following wc study the reductions for Lo^, L.i, Lr-,. which are solvable 
non-nilpotent Lie algebras. L3 contains an abelian ideal. L4 and L5 include a nilpotent 
ideal. 

2.3 Solvable algebra L3 = {T, P, i:>2, W} 

The algebra has the abelian ideal {T, P, VF}. Nonzero commutation relations are 

[A,T]=T, [P>2,P] = ^P. (2.66) 
The algebra has a decomposable structure 

L3 = {T,P,D2}©{iy}. (2.67) 
We note the canonical equation 

i^l^t + (1 + if2)^xx + (e + i7)|V'|'V' = (2.68) 
with the constants e = ^pl; /2 7^ 0, 7 and proceed to find the reduced ODEs. 

2.3.1 Subalgebra L3 1 = {T} 

A solution of VCNLS invariant under the algebra L3.1 must have the form 

i){x, t) = M{x) exp {iP{x)) . (2.69) 

M and P are found from the following reduced system 

eM^ - 2f2M'P' - M{P'^ + /2P") + M" = 0, (2.70a) 
7M^ + 2M'P' + /2M" + M(-/2P'^ + P") = 0. (2.70b) 

It is possible to arrange these equations as 

^ ^ ^•^^ - MP'^ + M" = 0, (2.71a) 



2 



1 + /. 

7 - ^4 ^ (^2p/y ^ Q_ (2.71b) 



2 



1+/. 

Similar to the preceding calculations, we were able to achieve decoupling for special 
cases of the parameters as follows. 
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(i) The case 7 7^ e/2. The function Y{x) defined as 



{M^P'y = J M* = Y' (2.72) 

1 + /2 

gives us the relations 

When we use these in fl2.71al) . we obtain a third-order equation for Y{x) 



Y'Y'" - h'- + ^^'±2ll\ U^Y- + ^^rf\ Y + CfY' = 0. (2.74) 



For a solution of this equation one can evaluate M and P from equations fl2.73p . 
(ii) The case 7 = e/2. 

M(a;) is a solution to the equation 

= ^ - eM^ (2.75) 

and P(x) is given by 

Pi^) = [ ^dx. (2.76) 



M2 



2.3.2 Subalgebra L3.2 = {P} 



Modulus and phase for the group-invariant solution corresponding to the algebra L3.2, 
which is in the form tp{x,t) = M{t) exp {iP(t))^ are found from the system 

eM^ - P' = 0, (2.77a) 
7M3 + M' = (2.77b) 



as 



M(t) = v/27t + Mi, (2.78a) 
P(t) =6(7^2 + Mit) + Pi. (2.78b) 



2.3.3 Subalgebra L3.3 = {T + eiVT} 

The solution in this case should have the form 



^(x,t) = M(x)exp i[eit + P{x)) , (2.79) 
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where M, P satisfy 



eM^ - 2f2M'P' - (ei + P'^ + f2P")M + M" = 0, (2.80a) 
+ 2M'P' + (-/sP'^ + P")M + /2M" = 0. (2.80b) 



In order to decouple these equations we can arrange them as 

'■-M^ - MP'"^ + M" = 0, (2.81a) 



g + 7/2 j^^3 , 
1 + /^ 



2 

(7 - ef2)M' + e,hM^ + (1 + fl){M''P')' = 0. (2.81b) 



Since /2 7^ 0, a first integral of fl2.81bl) can be obtained if 7 = e/2. Let Y{x) be 
defined as 

(M^P')' = — = Y'. (2.82) 



Then we have 

M^.-'-±llr, ^^-^^ (2.83) 
and thus fl2.81al) is transformed to an equation in terms of Y{x) 

Y'Y'" - -Y"^ - Ml±^y'3 _ J^1]-{Y + Cf = 0. (2.84) 



2 ei/2 1 + f, 



2 



2.3.4 Subalgebra L3.4 = {P + eiVT} 
In this case we have 



where 



We find that 



?/'(x,t) = M(t)exp i{tiX + P{t)) , (2.85) 



1 - eM^ + P' = 0, (2.86a) 
7M3 - /2M + M' = 0. (2.86b) 



^ -1/2 

7 ' 



M(t) = ( Ml exp(-2/2 + ^ J . (2-87a) 

^ln(Mi + ^exp(2/2t)) -t + Pi, 77^0 
P(t) = <^ ^ ^ (2.87b) 



2/2 Ml 



exp(2/2t) -t + Pi, 7 = 0. 
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2.3.5 Subalgebra L3.5 = {D + aW} 
We will look for the solution in the form 

1 



= -M(Oexp i{2a\nx + P{^) 



X 



X 

T 



(2.88) 



and the reduced system for M, P is 



eM^ + (2 + 6a f 2 - Aa^)M - 2(1 + Aaf2)iM' + (e' + 2(/2 - Aa)i)MP' 
-8f2^^M'P' - 4^2MP'2 + A^'^M" - Af2i^MP" = 0, (2.89a) 

+ (2/2 - 6a - 4a2/2)M + (8a - 2/2 - O^M' - 2(1 + Aaf2)^MP' 
+8^'^M'P' - 4/2^^MP'2 + Af2i^M" + Ai^MP" = 0. (2.89b) 

It is possible to rewrite these equations as 



(e + 7/2)M3 ^ 2(1 - 2a2)(l + /|)M - (2(1 + /2)e + 

+ (^2 _ 8^(^ ^ /2)e)MP' + 4(1 + f^)e{M" - MP'') = 0, 



e 



(7 - e/2)M^ - 6a(l + g)M' + (4a(l + ff^ ' y ) (M^)' 

-(2(1 + /2)^ + M')M2P' + 4(1 + /2)^2(^2p/y ^ Q 



(2.90a) 



(2.90b) 



so that they contain second order derivatives in terms of only M or P. But 
unfortunately we have not been able to proceed further as in the previous algebras. 

2.3.6 Subalgebra L3.6 = T + eiP + bW 

The modulus M and the phase P of the group-invariant solution 



V^(x,t) = M(Oexp i{bt + P{0) 



^ = X — tit 



(2.91) 



satisfy the system 



(e + 7/2)M3 _ 5m - ei/2M' + e^MP' - (1 + /|)MP'' 
+ (1 + fDM" = 0, 

(7 - ef2)M' + 6/2M2 - |(M2)' - ei/2M2p' 
+ (1 + /2^)(MV)' = 0. 

Arranging fl2.92bp according to the terms M'P' and we write 



2 ID'V 



{M'P 



1 + /. 



2(M2p') 



1 + /2' V 2 



^1 / ;\^2n/ 



(M^)' - bf2M' + (e/2 - 7)M^ 



(2.92a) 
(2.92b) 

(2.93) 
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We multiply this equality by exp (-[^^ define Y{^) as 



1 + 



First we have 



ei/2 



M2p' = exp(^0(^ + ^)- 
On the other hand, we need to solve for from 



(M'y-bf2M' + {ef2-j)M^)=Y'. 

(2.94) 

(2.95) 



ei/2 



(M^)' - 2e^bhM' + 2ei(e/2 - 7)(M2)2 = 2ei(l + Z^^) exp (-^ ^ (2.96) 

This equation is of Riccati type in and a special solution is needed for its 
integration. Furthermore, it can be linearized through the transformation = 
„ , \ — T^T at the cost of having its order increased by one: 

U" - 2e,bh U' + 4(7 - e/2)(l + fl) exp i)r U = 0. (2.97) 

Still, this equation does not lead to any immediate solution. Instead, it will be easier 
to handle fl2.96p by the choice 7 = e/2 

(exp ( - 2e,hM) M'J = 2e,{l + fi)exp (^^^(^^ - 26)^) Y' 
and if 6 = 2(1+/^) then we can integrate to find 



(2.98) 



\0 = 2ei(l + /I) exp {y{0 + c) • 



(2.99) 



A comparison of this result by (12.951) forces 



P' 



2ei(l + /|)- 



(2.100) 



Therefore we could end up with a decoupled equation for M from (I2.92ap 



M" 



l + fi 



4(1 + fiy 



M - eM\ 



(2.101) 



On the other hand, if (12.92bp is arranged with the condition 7 = e/2 as 

((1 + f2)M^P' - = ei^M^P' - 6/2M^ (2.102) 
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the choice b = 2(1+/^) ^ven makes it possible to write this equation in the simpler 
form 

U' = XU (2.103) 

where f/(0 = (1 + /^M^P' - f M^, A = The obvious solution f/(0 = 

Aoexp(A.^) with some constant Aq gives us the formula for P' 



Substitution of this result in fl2.92ap gives the decoupled equation for M as 

M" ^ ^Af + - eM3 + ^I^Pp) M-3. (2.105) 

This equation reduces to (12.1011) for Aq = 0. 







[P, D2] 


= M P 


[B,D2] , 


1 V B 



2.4 Solvable algebra L4 = {P, S, Da, W} 

This solvable non-nilpotent algebra is the extension of the nilpotent three dimensional 
Lie algebra {W,P,B}. We represent the action of D2 on this ideal by a matrix M 

/ 

M= 1/2 
\ -1/2 

We note that the algebra is not decomposable and the representative equation from 
Table 1 is 

#t + ^P.. + (e + tlM^^ + = 0, (2.106) 
where e = =f1, /i2 7^ and 7 are real constants. 

2.4.1 Subalgebra L4.1 = {P} 
Group-invariant solution is of the form 

^(x, t) = M{t) exp (zP(t)) (2. 107) 

and the reduced system is first order 

eM^ - P' = 0, (2.108a) 
M' + y M + 7M=^ = 0. (2.108b) 
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Pit) 



Solution of this system is elementary 

M{t) = (27tlnt + Mit)-i/2 

^ln(27lnt + Mi) + Pi, 77^0 

for /i2 = 1/2, whereas 

|^ln(Mi + ^ii-2'^^) + Pi, 7^0 

AMT^^'-^'^ + ^i' 7 = 

for h2 7^ 1/2. 



2.4.2 Subalgebra L4.2 = {B} 

The solution invariant under B has the form 

'ijj(x,t) = M(i)exp 

The reduced first order system is 

eM^ - P' = 0, 



which can be integrated to give 



2t 



a+2/i2 _ J_ j.\-i/2 



M(t)^{Mie+"'' --^ty 



Pit) 



^\n{M,-^r'^-)+P„ 7^0 



2/12 Ml 



7 = 0. 



2.4.3 Subalgebra L4.3 = {^^ + d^} 

The corresponding invariant solution is given by 



i^{x,t) = M{t)exp 

The reduced system becomes 

eM^ - P' = 0, 



4(l + eii) 



^Pit) 



t 2(l + eii) 
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M + 7M3 = 0. 



For the case /12 = 1/2 the system is integrated easily 

/ , \ -1/2 

M{t) = Kl + eit) (27 ln( + Ml n , (2.116a) 

r ^ln(27ln(^) + Mi)+Pi, 77^0 
P{t) = { ^ ^ ^ (2.116b) 

[wMTT^t) + Pi^ 7 = 0. 

U h2 7^ 1/2, then integration is possible in terms of the Gauss' hypergeometric 
function 2Fi{a, b, c; t): 



^1/2 



M{t) = (^Mit2'^^(l + eit) + Y^|^t(l + eit) 2^1(1 - 2/i2, 1, 2 - 2/12; -eit) 
P(t) = j eM^dt. (2.117) 

We note that for integer values of 2/i2 this solution will simplify to elementary 
functions. 

2.4.4 Subalgebra L4.4 = {D + aW} 

A group-invariant solution invariant under the subalgebra L4 4 will be of the form 



1 



^{x,t) = -M(Oexp ^i(2alnx + P(0) 



X 

T 



(2.118) 



Functions M, P will be solutions of the system 

eM^ + 2(1 - 2a^)M - 2^M' + (^^ _ 8a^)MP' - A^^MP'"^ 
+A^'^M" = 0, 

7M^ + (/i2^ - 6a)M2 + (4ae - y ) (M^)' - 2^M2p' 
+4^2(jvf2p'y = 0. 

For 7 = we define in (]2.119bl) the function V(^) as 



f2 

4^2(M2p')' - 2^M2p' = - 4a^)(M2)' + (6a - /i20M^ = 



From these relations we find 

M'^P 

and for /i2 = 1/4 



fl/2 
2 D' _ "5 



(2.119a) 
(2.119b) 

(2.120) 

(2.121) 
(2.122) 
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Thus if /i2 = 1/4 we can obtain from (KTTDf and fl2322|) that 



P(^) = I -aln^ + Pi. (2.123) 



This special form of -P(0 is readily seen to satisfy f l2.119bp whereas f|2.119ap is 
decoupled to determine M from 

M" = —M' -(—- — + -1-)M - -^M\ (2.124) 
2^ ^64 4^ 2^2^ 4^2 V ; 

2.5 Solvable algebra L5 = {P, fi, C2, W) 

L5 is another canonical extension of the nilpotent algebra {W,P,B} to a solvable 
non-nilpotent indecomposable four- dimensional algebra. The element C2 acts on the 
ideal {W^P^B} by the matrix M as 

W \ 

M = 



[W,C2] ^ 


\ ("^ 


[P, C2] 


= M P 




1 \ B 




Thus the last canonical equation under investigation will be 

#. + fe + i^l^f^H- ''%-;;'fg-'' ».-0. (2^125) 
with the constants e = =f1, /ii, /12 and 7. 
2.5.1 Subalgebra L5.1 = {P} 

The group-invariant solution ip{x,t) = M (t) exp (iP(t)) will be obtained by solving 
M, P from the reduced system 

' -M2 + -^-P' = 0, (2.126a) 



l + t2 ' l + t2 

Through the transformation M = W~'^^'^ in the Bernoulli type equation fl2.126bp we 
have 

(Vl + ^2 exp(-2/i2 arctant) w) = exp(-2/i2 arctant). (2.127) 

W can be expressed by a quadrature and hence 

, ^, , / exp (2/2,2 arctant) /, ^ /" exp(— 2/2,2 arctant) / 

Mit) ^ ^^^-^==^{m. + 2,J , (2.128a) 

/iiarctant+^lnfMi + 27/2iPi^^^^i^dt) +Pi, 77^0 
h arctant + ^ / ^-p^"^^^^"*) ^/t + Pi, 7 = 0. 
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We note that by the transformation t = tanz in the integrands these results can be 
expressed in terms of hypergeometric functions. 



2.5.2 Subalgebra L5.2 = {B} 
The solution has the form 

^Pix,t) = M{t)exp[t{^ + P{t)) 

with functions M, P determined from the system 



-M^ + - P' = 0, (2.129a) 



1 + ^2 '1 + ^2 

A transformation M = W~^^'^ in equation (]2.129bp leads to 

^ exp(— 2/j.2 arctant) = ■^--^=X== exp(— 2/i2 arctant). (2.130) 

We find M, P as 




t exp(2^2arctant) ^^^^^^j eM-^h^^nt) ^^^y"\^^^^^^ 



h, arctan t + ^ In (m^ + 27 / -P(-2j^^^'^^) dt) + P,, 7 7^ 

^ ^ ^ ^ ^ (2.131b) 

h, arctant + ^ / --p(-2fe^tan.) ^ ^ ^ q_ 

The transformation t = tan 2; is also applicable in the integrands leading to 
hypergeometric solutions. 

2.5.3 Subalgebra L5.3 = {C2 + aW} 
Group-invariant solution must have the form 

ij{x,t) = M(Oexp [z(a arctant + + P(0)1, ^ = J—:, (2.132) 

L V 4(1 + V) /J Vl + 1 

Substitution of this solution into the original equation ends up with the system 

eM^ + {hi-a-^)M - MP'^ + M" = 0, (2.133a) 
7M^ + /isM^ + (M^P')' = 0. (2.133b) 
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In the case 7 = 0, /i2 7^ decoupling of these equations is possible 



Y'Y'" - -Y"^ - ^Y'^ + {2hi -2a- ^)Y'^ - 2hl{Y + Cf = 0, (2.134) 
2 112 2 

M=(-lr)V2, p' = -h2^^. (2.135) 

n,2 i 

If 7 = /12 = then we have a second order nonlinear ODE for M 

M" = C^M-^ + (^ + a - hi)M - eM^ (2.136a) 

p' = CM-"^. (2.136b) 



3 Analysis of the reduced equations 

In Table [3] we refer to the numbers of the reduced system of equations of first order 
besides the second and third order equations obtained through the decoupling task 
for some special values of the arbitrary parameters. Among the third order equations 
we only included those which are polynomials in the derivative Y' , e.g. f l2.18p is not 
in the list. We have expressed the solutions of first order equations in the preceding 
Section. This part of the work will be devoted to the study of solutions of the second 
and third order equations. 



Table 3: Equations under study 



Order 


Equation Number 






1 


fl2.46p. (12.631). 02.77P. fl2.86p. 
02.1151). fl2.126p. 02.129P 


fj2.1()8p 


, mm, 


2 


(I2.2ip. (12.341). fl2.42p. fl2.75p. 


fl2.105p 


. 02.124P. fl2.136p 


3 


02.13P. (12.281). 02.39P. fl2.53p. 


fl2.60p. 


fl2.84p. fl2.134p 



3.1 Third order equations 

None of the seven third order equations summarized in Table [3] passes the invariant 
Painleve test for PDEs. Since fl2.53p and f l2.84p do not contain the independent 
variable, we can directly lower their order by one if we set Y' = W{Y). We obtain 
the following second order equations: 
(i) For equation flZ^ with W = 

1 ■ . 2e 2a 2hl{Y + Cf , , 

W = ^ \ \ —. (3.1) 

2W h2 W ^ ^ 
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(ii) For equation fl2.84p . 

2W ^ (1/2 1 + ^ ' 

We seeked a first integral to these second order equations in the form 

A{Y, W)W'^ + B{Y, W)W'^ + F{Y, W)W' + G{Y, W) = I (3.3) 

for some functions A, B, F, G and a constant /, but saw that it is not possible. 

For all the third order equations satisfied by Y = Y{^) (including f l2.53p and 
f l2.84p ) we suggested a first integral of the form 

A{^, Y, Y') Y"' + i?(e, Y, Y') Y" + F{^, Y, Y') = I (3.4) 

with some functions A, B, F and a constant I. A first integral of this type exists only 
possible for fl2.39p for the special values of the constants c = 0, hi = {5 + 81hl)/36. 
The first integral has the form 

+ -'^Y'' + + ^^^^ Y" + ^ YY' (3.5) 

3a; h2 9x^ 

,l2hlC I , Ahl .^^ ^ 
+ 1 jj^)Y' + ^{Y + Cf = {). 

X'^ x'l-^ x^ 

We made an unsuccessful attempt to find a further first integral 

A[x, Y) + B{x, Y) Y'^ + F{x, Y) Y' + G{x, Y) = K (3.6) 
where K is another constant. 



3.2 Second order equations 

Among the second order equations successfully decoupled from the reduced systems, 
fl2.42p passes the P-test for hi = 5/36 and so does fl2.75p without any condition on 
the parameters. For the other five equations which do not pass the P-test, we looked 
for a first integral, assuming M = M(^) satisfies an equation of type 

A(^, M)M'3 + B{^, M)M'^ + F(^, M)M' + G{^, M) = /, (3.7) 

which ended up without any success. 

Before proceeding to the solutions of second order equations passing the P-test, 
we note that equation f l2.105p does not contain the independent variable if we choose 
Ao = 0, which means a reduction in order. Indeed, if we set M' = aiW{M) with 
ai = Tr%5 an Abel equation of the second kind is obtained 

Vr(Vr - 1) = - — M^, (3.8) 
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where 02 = ^dlf'i) ■ For n = ^ — 3 and 62 = T^, w = w{z) a. transformation in 

the parametric form 



n+2 1 n + 2 , , 71 -\- 2, ^ , , 

M = 2~w, W = z~ (zw' -\ w) (3.9) 

n + 3 2 

converts this equation with A = ——^7x573^ to an equation of Emden- Fowler type E] 

w" = Az^'w^ (3.10) 

which drove the final nail in the coffin. 

We close this section with the analysis of equations passing the P-test. Painleve 
and his successors classified second order differential equations that have at most 
pole-type singularities in all their solutions and determined such 50 equivalence classes 
together with their representative equations. For details the interested reader is 
referred to [5]. Since we have two second order equations passing the P-test, we 
are going to try to find the equivalence class to which they may belong. 

Transformation of fl2.75p to the equations numbered PXVIII and PXXXIII in the 
Painleve classification of second order nonlinear ODEs, their first integrals and hence 
solutions in terms of elementary and elliptic functions in various cases were done in 
|10] . Since a simple substitution and a careful account of the different cases depending 
on the constants will suffice to find the results for ( 12.75 p . we do not reproduce them 
here and refer the interested reader to that work. 

There remains the treatment of equation fl2.42p . If we make a change of the 
dependent variable as M = a/ H{x), H{x) > we have 



(3.12) 



We apply a further transformation H{x) = X{x)W{ri{x)) and find 

2eA . 2C2 1 

W H W . 

fj'^x X^jf 

We will determine A, i] such that this equation is of Painleve-type. 
(i) The case C 7^ 0. If A, 77 and other constants are chosen as 

V = r,ox^'\ A = Aoxl/^ c<0, r/o = -(y)'^', ^^i = ^ (3-13) 
an equation quite similar to PXXXIV is obtained. 

W = + AaW^ - 'qW + , (3.14) 
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where a = —^2^2) ^ — ^^'^ -^o is arbitrary. By a final 
transformation 

2aW = V + V'^ + ^ (3.15) 
we see tliat V satisfies the equation 

V = 2V^ + T]V + k, k = -]^±Aa5i. (3.16) 

which is the second Painleve transcendent so that we can express V = Pii{r]o x'^^^). 
Since W is complex-valued Aq has to be chosen so that the product XW is real, 
(ii) The case C = 0. If we choose 

V = Vox'/\ A = Aoa;V3, r/o=(^)'^', ^o = -e(^)'^', h = ^ (3.17) 
we arrive at the equation PXX: 

W = + + 27]W. (3.18) 
2W 

Setting = W leads to Pu agam 

i) = 2U^ + 7]U. 
We can explicitly give the solution as 

^ = X'J^x^/'Pn{r]ox^/') exp{i{ct + Pq))- (3.19) 

4 Solutions by truncation 

In order to investigate Painleve property for (11.11) in [1] we wrote the equation with 
its complex conjugate as the system 



iut + f{x, t)uxx + g{x., t)u^v + h{x, t)u = 0, 
—ivt + p(x, t)vxx + 9(3;, t)uv^ + r(x, t)v = 



2 , . (4.1) 



and expanded m, f as 

00 00 
u{x,t) = J2 Uj{x,t)<^>''+^{x,t), v{x,t) = J2 Vj{x,t)<l>'^+^{x,t). (4.2) 

j=0 j=0 

We obtained the conditions on f,g,h so that all solutions to VCNLS are in this 
form. Coefficients for canonical equations of 4-dimensional subalgebras do not satisfy 
the compatibility conditions of the P-test, therefore they do not have the Painleve 
property. Since the conditions obtained from the P-test are equivalent to those for 
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having a Lax pair, they are not integrable OH]. In this case if the series (14. 2 p is 
truncated at an order j = N and plugged in the equation, a system of equations for 
Wj, J < and $ has to be satisfied. An exact solution is obtained once this system 
can be solved in a consistent way. 

For the Painleve test to be successful it is required that resonance coefficients uj 
corresponding to the resonance indices j = —1,0,3,4 are arbitrary. This is true if 
the compatibility conditions at resonance levels hold. As we already mentioned, this 
is not the case for our canonical equations, and we first checked whether resonance 
equations are satisfied at all for some special form of u/s and $. The results were 
not so promising, since either no condition for $ or conditions being equivalent to 
the integrable case can arise. When we were lucky to obtain a specific form for $, 
conditions other than resonance levels did not hold. Therefore, we could not obtain an 
exact solution and a Backlund transformation by the truncation approach. However, 
when we applied the method as it was done in [TT], we were able to obtain nontrivial 
exact solutions. 

As the first step of the Painleve test the leading orders a and /3 are determined 
by substitution of u ~ mo$" and v ~ fo$^ in (14. ip . Balancing the terms of smallest 
order requires that 

a + /3 = -2 (4.3) 

and 

uovo = -a{a - 1)^$^. = -/3(/3 - 1)^$2 (4.4) 

9 q 

hold. Since the leading orders should be negative integers for the equation to have the 
Painleve property, (14. 3 p implies that a = —1 and /3 = — 1. Since we are interested in a 
case in which the equations do not have the P-property, we weaken the condition that 
a and /3 are integers and determine the leading orders by solving the eqs. (14. 3 p and 
(14. 4 p simultaneously. This will indeed lead to finding exact solutions by truncation 
approach. 

We applied this approach to the canonical equations of algebras Li, L3, Z/4 
successfully. Overdetermined system of equations for L2 and L5 algebras are not 
compatible and the method fails to apply. 



4.1 Truncation method for the algebra Li 

The coefficients for the algebra Li are f = 1, g = {e + h = {hi + ih2)^. We 

apply the truncation method to the slightly more general coefficients 

/=!, g = {e + t^)^, h= {hi + th2)—^, a,beR. (4.5) 

If a 7^ 1, 6 = 2, the equation is invariant under the 3-dimensional solvable algebra 
with a basis 
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The algebra is extended for a = 1. 

When we solve (14.31) and fl4.4p together, for 7 7^ we have 



a 



-1 - i5, 13 



-l + i6: 



6 



-3e ± + 9 

2^7 



and (14. 4p simplifies as 



uqVo = a;"<l>^.. 

7 



(4.6) 



(4.7) 



We truncate the Painleve expansion at the first order (j = 0) and suggest that solution 
has the form 



-1-iS 



v{x,t) = fo(x, t)$(x, t)' 



-l+iS 



Putting these expressions in (14.11) . the terms <|)-3±i<5^ ^-2±t(5^ ^-i±j<5 ^jjj appear. We 
choose the coefficients of these terms equal to zero to obtain a system of three 
equations for uq, Vq and $ each of which consists of two equations. The condition 
obtained at the level ^-^i*"? jg (|4.7p itself. Coefficients of ^-^i*-? vanish if 



2$. Mo 2$ 



0, 



+ 



2$^ vo 2$^. 
Taking the difference of these two equations and using (14.70 lead to 

0o(t)xi-t +0i(t), a^3 
(poit) In |x| + 01 (t), a = 3 



(4.9a) 
(4.9b) 



$(x,t) 



(4.10) 



with arbitrary functions 0o(^) and 0i(t). We can also solve Uq, Vq by using this form 
of $ in (USD. 



Uo{x,t) = Ci{t) x"/^ exp 
Vo{x,t) = C2{t) x"/^ exp 



2$. 
ax 



■dx 



(4.11a) 
(4.11b) 



Here Ci{t) and C2{t) are the arbitrary constants of integration and the expression in 
the exponential function can be evaluated as 



2$. 



-dx 



IX 

4 



f{\n\x\ 



^(iox^ _l_ <i>i In k'l N 
n</,o 8 00 4 ^ 



4{l-f) ' <^0 2(1-4) 



a = —3 



(4.12) 
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Functions u and v are treated independently in applying the P-test. But in this case, 
in order that the system (14. ip corresponds to (II. ip . we must have v* = u. This 
requires that C2{t)* = Ci{t). As a conclusion of this condition, we get 

\Ci{t)\'={ ^ (4.13) 

using the special singularity (I4.10p in equation (14. 7p of order <|)-3±i(5_ 

Again the constants 5 and 7 must satisfy ^ < 0. This means that we have to 
choose the negative sign for the formula of 6 in (14. 6 p 
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We solve the equations obtained by choosing the coefficients of terms $^1='="^ zero 
in various cases depending on the constants a, b. 
(1.) The case a = 3. 

We solve the equations at order $~i='="^ and find b = 2, hi = 1/4, h2 = 0. As a 
result we have 

Uo{x,t) = c^/x, <l>(x, t) = fco In |x| + /ci, (4-15) 



where ko, ki are real arbitrary constants and c G C with |c| = ko^/ — ^. We write the 
solution explicitly as 



u{x, t) = - — — ^ — — exp ( — z5 ln(fco In |x| + fci)) . (4.16) 

(2.) The case a = —3. 

For the constants we have b = 2, h"^ = 3 + 4hi. uq and $ are found to be 

Mo(x,t) = cx"5+*^, ^{x,t) = kox'^ -2h2kot + ki, (4.17) 



X 

"(^' ^) = -UTT, — 9 \, , ^ , , ^ exp ( i In - — ^, ■ ^ . ) . (4.18) 



ko, ki are arbitrary real constants and c G C with the modulus |c| = 2kQ^/ — ^. The 
solution is given by 

(A;ox2 - 2/12M + ^1) ^^^'^ V'"' (^oa:^ - '^hk'ot + ^i)^ 7 ' 

(3.) The case a 7^ =f3. 

The system of equations look quite complicated in this case. However, 
compatibility implies that one of them can be expressed in a relatively simple form 

\Gi G2 a — 6 (po I a — 3 00 
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In order to complete the analysis, we need to find for different values of a, b the 
functions Ci, 00? 0i that will satisfy fl4.19p and the accompanying condition 

a(a-6) fC^ , 9 2^+, 

+ f4^i - -T^'^) (4-20) 



2(a-3)2 02 2(a-3) 
18 



(a + 3)(a-3)2 



^ 00 00 



We proceed by considering powers of x in ( I4.19p . 
(3.i) The case = 6 = — 1 + | + We require a = 3. This case is not possible since 
we had been able to find an exact solution in (1.) for 6 = 2. 

(3.ii) The case 07^6=— 1 + | + 6. Since we must have a = 3, this case corresponds 
to (1.). 

(3.iii) The case = 67^— 1 + | + 6. Equation f l4.19p requires a = or 0i = 0. 

(A) The case a 7^ 0. We find 0i(t) = fci as a constant. Since the coefficient of in 

fl4.2Up must vanish, we find a = 6. (14.201) and fl4.19p simplifies to 

'Ci C-i 



— + — 

C\ C2 (po 



2/^2 + ^ + ^ + 3^ = 0. (4.22) 



Coefficient of in (I4.2ip gives 

Mt) = ht + h (4.23) 

with constants ko, and integration of (I4.22p results in 

exp(2/i2t) Ci C2 00 = s, s = const. (4.24) 

Remembering that C2 = C* and considering |Cip = — ^0o from fl4.13p . we see 
that fl4.23p and f l4.24p are satisfied only if h2 = k^ = 0. When we set Ci{t) = 

koJ-—exp (iK{t)), (I4.2ip implies that K{t) = hit + k2. Eventually we have 



uo{x,t) = ko^l X exp ^i(/iit + ^2) j , ^{x,t) = kox ^ + ki (4-25) 

and 



I 36 x'^ / \ 

u{x,t) = koJ — exp (i(hit -5ln{koX~^ + ki) + ^2) ). (4.26) 

y 7 + kix V / 
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(B) The case a = 0. Through a similar analysis we have (f)o{t) = ko, 0i(t) = kit + ^2, 
Ci{t) = ko.r^ exp {iK{t)), K{t) = {hi - + k^ together with the condition 



h2 = 0. Therefore we get 



I 36 \ ( k^ ki \i 

Uo{x,t) = koi exp ii{hi--^)t--—x+k3) , ^{x,t) = kox+kit+k2 (4.27) 

V 7 L V 4fcs 2kn /J 



7 ^ IV ' 4k',' 2ko 
and the explicit solution 



.35 ^2 ^ 

^(^'^) = 1 "X /, , exp [z('(/ii--L)t--^x-51n(/coX + M + fc2) + A:3)l (4.28) 

knx + k-it + ko I V 4A;n 2kn /J 



is obtained. 

(3.iv) The case 67^0 = — 1 + | + There is no solution in the truncated expansion 
form. 

(3.v) The case 6 0, -1 + f + 6 7^ 0, 6 ^ -1 + f + 6. All the coefficients in flCTjl 
must vanish: 



h2 = 0, (4.29a) 

CiC2<Pp = s, (4.29b) 
a<pi = 0. (4.29c) 

If we remember that C2 = and that = _ |)202 ^^^^ (KT3\f . K2SUf 

becomes 

- -(1 - f = s. (4.30) 

Analysis should be carried out for three different cases according to the constant a 
under the conditions fl4.29cl) and fl4.30p . 

(A) The case a = {b ^ {0,1}). flTOj) and (^M> together require hi = hi = 0, 
(poit) = ko, 01 (t) = kit + k2, Ci{t) = cexp(-§t), |cp = -^^O' c G C. Functions 
Mo, $ generate 

Mo(a;, t) = cexp ( - ■i(— — X + — ijt) ), ^(x,t) = kox + kit + k2 (4-31) 
V ^2ko 4kQ / 

and therefore the exact solution is given by 

u{x,t) = - ————exp -i(—x+-Y2t + 5\n{koX + kit + k2)). (4.32) 

kox + kit + k2 L \2A;o 4A;o /-I 

(B) The case a ^ {0, 1}. We find from and dOO]) that 0o(t) = ^o, 0i(t) = ^1, 
^1 = ^2 = 0, a = 6, Ci(it) = c and |cp = — ^^o- We conclude that 

Mo(a;,t) = cx, ^{x,t) = kox^^ + ki (4.33) 
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and 

cx^ f \ 

u{x, t) = — exp - i5 ln(A;oa;~^ + ki)]- (4.34) 

kq + kix V / 

(C) The case a = 1 (6 7^ 2/3). This last situation is going to give us the exact solution 
for the canonical equation of the algebra Li, namely the case a = 1 and 6 = 2. 

Equation (14.19P requires /i2 = and (pi = ki. f l4.20p is satisfied if and only if 6 = 2 
and hi = 5/36. fOOll and f lCTD take the forms 



2 

^^' + 7(t^-^^)^' = 0' (4-35) 



Ci _ C2\ . 2 , 3 /0O _ 5 00 \ ^4 

Ci Cj'"' 4^00 2 



Ci C2 00 



x^ = 0. (4.36) 



We arrive at the functions 0o(t) = {kot+k2) ^ , Ci{t) = Ci0o(t) and C2{t) = C20o(t). 

-AS 
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It is obvious that C2 = and the condition fl4.13p requires |cip = ci G C. We 



sum up the results so far: 



= WTW l^ 4(M + fc2) >' ' ^'-'^^ 

^x, t) = {kot + A;2)"^/^ x^/^ + ki (4.38) 



and hence 



Remark. Since the canonical equation of Li is invariant under the action of 
the group of transformations S'L(2,M), which is the composed action of translation 
generated by T, scaling Di and the conformal transformation of Ci, the solution 
f l4.39p is transformed into a new solution of the canonical equation under this action. 
Owing to the invariance property under Ci, this transformed solution has a finite 
time singularity and therefore was fruitful to study blow-up profiles. It is exactly 
this blow-up character that was used in [12] to establish the existence of singular 
behaviours of solutions in the sense of Lp and L^o norms and in the distributional 
sense as well. 

4.2 Truncation method for the algebra L3 

Coefficient functions for the algebra L3 are 

/ = l + z/2, g = {e + i^), h = 0. (4.40) 

In fact, this constant coefficient case is included in [11] but we could not deduce our 
results from theirs. Differing from the previous algebra, / contains imaginary part 
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and there will be a slight modification in the above construction. If we solve (14. 3p 
and (14. 4p together we find the leading orders to be 

-3(e + 7/2) ± v/9(e + 7/2)2 + 8(7-6/2)2 



a = —1 — i5, (3 = —1 + iS; 5 



for 7 7^ e/2 and (14.41) is equivalent to the condition 



2(7 - e/2) 



(4.41) 



novo = -^^^^^6^1. (4.42) 
7 - e/2 

We truncate the Painleve expansion at the first term and propose a solution of the 
form (14.81) . We try to determine Uq, Vq and $ by choosing the coefficients of the terms 
^-3±i5^ ^-2±iS^ ^-i±iS Tf^iiiQi^ appear when these ansatze for u and v are put in (14. ip . 
Equations for order $-^±1(5 equivalent to (14.421) . Terms $-^±1(5 (disappear if 

+ 2(1 + z/2)^ + (1 + z/2)|p = 0, (4.43a) 

Uo 

+ 2(1 - ^f2f-^ + (1 - ^/2)^ = 0. (4.43b) 
These equations can be arranged as 



2(1 + /I) Mo 2$, 

£i^_L + Jli^ + ^£ = 0. (4.44b) 

2(1 + /I) t;o 2$, ^ ^ 

Adding equations (14.430 and integrating with respect to x and doing the same for 
equations (I4.44p result in 

uoVo( — ) $x- = i^i(t) (4.45a) 

exp j ^dxy^v^^, = F2{t) (4.45b) 

with some arbitrary functions Fi and F2. Moreover, the fact that Vq = Uq brings 
f02|) to the form 

2^_3(1^^^^2^ (4.46) 
7 -e/2 

When we put this expression in the absolute form of (I4.45ap . we see that $(x,t) = 
00 (^) X + 01 (t). This leads to the result that the exponential term in fl4.45bp cannot 
depend on x. Therefore we have = 0, that is, $(x, t) = kox + ki with the constants 
ko, ki. Equations for the order ^^^'^'^^ become quite simple and we solve them to find 
uo{x,t) = c G C, |c|2 = — ^^^^/fj S fco- is necessary to choose the negative sign for 
the square root in the formula 6 of (14. 41 p . The corresponding exact solution will be 

^(2^) t) = 1 —r- exp ( -i6 \n(ko x + ki)] . (4.47) 

k^x + ki \ / 
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4.3 Truncation method for the algebra L4 

We repeat the arguments which worked for algebra Li for the potential h{x, t) = ^ 
of algebra L4 and find 

2 

«(^'^) = :rT17T- exp -SH-^, + h))] (4.48) 



X + kokit 



At kot 



if /12 = 1/2. Here the constant c G C must satisfy |c| = in addition we have 

^ _ -36-^/872+9 
~ 27 

References 

C. Ozemir and F. Giingor, On integrability of variable coefficient nonlinear 
Schrodinger equations. arXiv:1004.0852 v3 [nlin.SI]. 

Gagnon, L., Winternitz, P., 1993, Symmetry classes of variable coefficient 
Schrodinger equations, J. Phys. A:Math.Gen., (26) 7061-7076. 

Al Khawaja, U., 2010, A comparative analysis of Painleve, Lax pair, and 
similarity transformation methods in obtaining the integrability conditions of 
nonlinear Schrodinger equations. Journal of Math. Phys., 51, 053506. 

Brugarino, T., Sciacca, M., 2010, Integrability of a nonlinear Schrodinger 
equation. Journal of Math. Phys., (51), 093503. 

Olver, P.J., Applications of Lie Groups to Differential Equations. 
Springer- Verlag, NewYork, 1986. 

Ovsiannikov, L.V., 1983, Group analysis of differential equations. Academic 
Press, New York, 1982. 

P. Winternitz, 1992, Lie Groups and Solutions of Nonlinear Partial Differential 
Equations, in "Integrable Systems, Quantum Fields, and Quantum Field 
Theories, L. A. Ibort and M. A. Rodriguez, Eds." , Netherlands, Kluwer Academic 
Publishers. 

Polyanin, A.D., Zaitsev, V.F., 1995, Handbook of exact solutions for ordinary 
differential equations, CRC Press, USA, p44. 

Ince, E.L., Ordinary Differential Equations, Dover, New York. 

Ozemir, C., Giingor, F. Group-invariant solutions of the (2+l)-dimensional cubic 
Schrddinger equation, J. of Phys. A (39), 2973-2993. 



33 



[11] Yomba, E., Kofane, T.C, 1996, On exact solutions of the generalized modified 
Ginzburg-Landau equation using the Weiss- Tabor-Carnevale method, Physica 
Scmpta, (54), 576-580. 

[12] F. Giingor, M. Hasanov, C. Ozemir, Unimodular invariant variable coefficient 
nonlinear Schrodinger equation and blow-up of its solutions. larXiv:1101.2307l 
[math.AP]. 



34 



